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Abstract. We studied the colliding galaxy NGC7714 with two XMM-Newton observations, six months apart. The galaxy 
contains two bright X-ray sources: we show that they have different physical nature. The ofF-nuclear source is an accreting 
compact object, one of the brightest ultraluminous X-ray sources (ULXs) found to date. It showed spectral and luminosity 
changes between the two observations, from a low/soft to a high/hard state; in the high state, it reached Lx » 6 x 10'"' erg s"'. 
Its lightcurve in the high state suggests variability on a » 2 hr timescale. Its peculiar location, where the tidal bridge between 
NGC7714 and NGC 7715 joins the outer stellar ring of NGC7714, makes it an interesting example of the connection between 
gas flows in colliding galaxies and ULX formation. The nuclear source (L^ w 10"" erg s"' ) coincides with a starburst region, and 
is the combination of thin thermal plasma emission and a point-source contribution (with a power-law spectrum). Variability in 
the power-law component between the two observations hints at the presence of a single, bright point source (L^ S 3 x 10* erg 
s"'): either a hidden AGN or another ULX. 
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1. Introduction 

XMM-Newton and Chandra studies of colliding or merging 
gas-rich galaxies at various stages of the evolutionary sequence 
(Toomre 1977) have revealed significant contribution to the X- 
ray emission both from diffuse hot gas, associated to starburst 
processes, and from accreting point sources, generally associ- 
ated to a young stellar population. See, for example: the Mice 
(Read 2003); the Antennae (Zezas et al. 2002; Fabbiano et 
al. 2003a); M82 (Grifliths et al. 2000); the Cartwheel (Gao 
et al. 2003). A peculiar feature in many of these systems 
is the presence of accreting X-ray sources brighter than the 
Eddington limit for a stellar-mass black hole (BH) (LEdd ~ 10^^ 
erg s '); they are commonly known as ultraluminous X-ray 
sources (ULXs). The ages and masses of the compact objects in 
ULXs, the nature of the donor stars, and the geometry of emis- 
sion are still unclear, and hotly debated. It is also unclear pre- 
cisely why ULXs are preferentially found in interacting galax- 
ies (Swartz et al. 2003), and what this can reveal about their 
mechanism of formation. 

The interacting system Arp 284 (Arp 1966) is an excep- 
tional example of a recent (~ 100-200 Myr ago), direct im- 
pact (Struck & Smith 2003). It consists of the nuclear starburst 
galaxy NGC 7714 (classified as SB(s)b pec') and its fainter. 
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currently inactive companion NGC 7715 (Im pec). NGC 77 14 
is located at a redshift distance of 37 .3 Mpc (Huchra et al. 1 999, 
for Hq = 75 km s"' Mpc"'). It has a prominent stellar ring, 
three tidal arms/tails, and is connected to NGC 7715 by a gas 
and stellar bridge (Arp 1966). Its low inclination (viewing an- 
gle of 45°) and low foreground absorption (hh = 4.9 x 10^" 
cm"^; Dickey & Lockman 1990) make it a good target for X- 
ray studies. 

ROSAT/HSl observations of NGC 7714 (Papaderos & 
Fricke 1998) have revealed two strong (Lx > 10**" erg s"') X- 
ray sources, separated by ^ 22". The brighter one coincides 
with the starburst nucleus; the fainter one is located approx- 
imately where the gas/stellar bridge joins a collisional stellar 
ring, but does not correlate with any bright counterpart at other 
wavelengths, nor is it located in a starburst region. It was sug- 
gested (Papaderos & Fricke 1998) that the off-nuclear source 
might be a compact region of hot shocked gas. This could be 
due either to the collision of a fast starburst-driven outflow with 
the colder, denser gas in the galactic bridge; or, to the infall of 
high-velocity H I clouds along the bridge onto the outer H I disk 
of NGC 7714 (a somewhat similar situation to accretion-disk 
hot spots in X-ray binaries). However, the limited wavelength 
coverage and resolution of ROSAT did not allow detailed indi- 
vidual analyses of the two sources. 

Here we present some preliminary results of our XMM- 
Newton study of the system. We argue that the off-nuclear 
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Fig. 1. Top panel: true-color image of the NGC 7714 field ob- 
tained by coadding the MOS and pn observations of 2002 Jun 7 
and Dec 8. The image was smoothed with a 3 x 3 pixel boxcar. 
Red: 0.2-1 keV; green: 1-2 keV; blue: 2-12 keV. Size of the 
image: 7!8 x 6.'2; North is up. East is to the left. Bottom panel: 
on the same scale, the positions of the three brightest X-ray 
sources in the field are overplotted (red circles) onto a Digitized 
Sky Survey B image; the positions of a few other, fainter X-ray 
sources detected with XMM-NewtonfEPlC are also overplotted 
(green circles and ellipse). 

source is an accreting point source (in fact, one of the most 
luminous ULXs ever detected), and we discuss the spectral and 
timing properties of the two sources. 

2. Observations and data analysis 

NGC 77 14 was observed with all instruments on-boardXMM- 
Newton on 2002 June 7 and 2002 December 8 (Table 1); we 
used the medium filter, full-frame mode for the EPIC detectors. 
We processed the Observation Data Files with standard tasks 
in Version 5.4 of the Science Analysis System (SAS). After in- 
specting the background fluxes, we rejected the initial a; 1 ks 
of the first exposure, which was affected by a background flare: 



X-ray source in the colliding galaxy NGC 77 14 




Fig. 2. Position of the ULX and of the X-ray nucleus (red cir- 
cles with radius = 3") overplotted on an HTS/WFPC2 image in 
the V band (f555w filter). The ULX is located where the tidal 
bridge meets the coUisional steUar ring. North is up. East is to 
the left. 

we retained a live-time good-time-interval of 15.9 ks (for the 
MOSs) and 14.3 ks (for the pn). From the second observation, 
we used a live-time good-time-interval of 15.8 ks (MOSs) and 
12.9 ks (pn). We filtered the event files, selecting only the best- 
calibrated events (pattern < 12 for the MOSs, pattern < 4 for 
pn), and rejecting flagged events. We checked and corrected the 
astrometry of the various images so that the position of the nu- 
clear source coincided with the NED position (Clements 1983): 
R.A. (2000) = 23'' 36"' 14M, Dec. (2000) = -h02°09' 18'.'6. 

Firstly, we extracted source and background spectra for the 
nucleus and the bright off- nuclear source in NGC 7714, for 
each detector in each of the two exposures. The ofi'-nuclear 
source is located at R.A. (2000) = 23* 36'" 15'.6, Dec. (2000) = 
-h02° 09' 23'.'5, ie, ^ 22" east of the nucleus (Figs. 1,2). We es- 
timate an error % +3", mostly due to the astrometry uncertain- 
ties; however, the error in the relative positions of the nuclear 
and off-nuclear sources is only » 1". To reduce the contami- 
nation from the bright nearby nucleus, we used as source ex- 
traction region a circle of radius 11". We used as background 
region the union of three circles of 11" radius, centred at a 
distance of 22" from the nucleus, on the opposite side of the 
source, so that it would contain a similar contribution from the 
nuclear source. For the nucleus, we used an extraction circle 
of radius 15". The background region was chosen as the union 
of two circles of radius 15" centred at 22" from the nearby 
off'-nuclear source. We tested the efffect of different choices of 
background regions, until we were satisfied that the reciprocal 
contaminations of the two sources were properly subtracted. 

We built response functions for each source with the SAS 
task rmfgen, and auxiliary response functions with arfgen. 
Given the small extraction radii used for the two overlap- 
ping sources, we corrected for the ffux that is outside the 
extraction regions by setting "extendedsource=false", "mod- 
elee=true" in the arfgen parameters. We then fitted the 
background-subtracted spectra with standard models in XSPEC 
v.11.3 (Arnaud 1996); owing to the uncertainties in the EPIC 
responses at low energies, we used only the 0.3-12 keV range. 
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Table 1. XMM-Newton/EPlC observation log for NGC 7714. 3. Main results 



Date Obs ID Start time Stop time GTI 

(ks) 



2002Jun07 0112521301 pn: 10:36:39 15:16:39 15.8 

MOS: 10:03:23 15:21:23 16.1 

2002 Dec 08 0112522601 pn: 12:40:16 16:38:19 14.3 

MOS: 12:18:15 16:43:15 15.9 



3.1. Starburst Nucleus 

The nucleus of NGC 7714 is characterised by an active star- 
burst that started ~ a few xlO^ yr ago, with an average star- 
formation rate x IMq yr"'; the rate has probably been higher 
by a factor of a few during the most recent burst which started 
X 5 Myr ago (Langon et al. 2001). Its optical, UV and IR prop- 
erties are similar to those of the "prototypical" starburst galaxy 
M 82. In analogy with other starburst galaxies, we expect the 
nuclear X-ray emission to be a combination of diffuse hot gas, 
high-mass X-ray binaries and supernova remnants — including 
perhaps some contribution from the Type Ib/c SN 1999dn, lo- 
cated x 13" south-east of the nucleus (Qiu et al. 1999; Deng et 
al. 2000). 

The limited spatial resolution of XMM-Newton does not 
allow us to resolve individual sources in the nuclear region. 
However, we checked whether the observed profile is con- 
sistent with the point-spread-function of a point source, us- 
ing standard tasks in IRAF. We compared the full-width half- 
maximum (FWHM) of the radial profile for the nucleus, the 
off-nuclear source and a nearby quasar Using EPIC/MOS im- 
ages from the two observations, we obtain that the last two 
sources have FWHM » 6", while the nucleus has FWHM ^ 8". 
This suggests that the X-ray emission from the nuclear region 
is indeed extended, in agreement with the ROSATjW'Rl results 
(Papaderos & Fricke 1998). No difference in the radial extent 
of the three sources is found in the EPIC/pn images; however, 
this is due to the lower spatial resolution and larger pixel size 
of the pn. We also examined the radial profiles separately for 
"soft" (0.2-1.5 keV) and "hard" (1.5-12 keV) MOS images. 
We obtain that, in the hard band, the FWHM is the same for all 
three objects, consistent with point sources. The larger radial 
extent of the nuclear source is found only in the soft band, sug- 
gesting a different origin for the two components (see Section 
4.3). Further investigation of the MOS data suggests that the 
larger radial profile of the nuclear source in the soft band is due 
to faint emission a few arcsec south-west of the central posi- 
tion, in excess of what would be expected from a point-source 
point-spread-function. No excess emission is detected east of 
the central position, suggesting that SN 1999dn does not con- 
tribute significantly. 

We fitted the nuclear X-ray spectrum with an absorbed 
two-temperature variable-abundance thermal plasma model 
(vmekal), plus an absorbed power-law component, which ac- 
counts for the contribution from accreting compact objects. We 
find that we cannot fit the EPIC spectra from 2002 June and 
December simultaneously with the same set of parameters. The 
two spectra are identical at energies % 1 keV, but the flux at 
higher energies was higher in December Therefore, we fixed 
the two thermal-plasma components (which we do not expect 
to vary on a 6-months' timescale) but allowed the power-law 
contribution to vary between the two observations. This results 
in a good fit (xl - 0.77, see Table 2 and Fig. 3). We do not ob- 
tain statistically-significant improvements to the fit by allowing 
different normalisations of the thermal plasma components in 
the two observations. The total emitted luminosity in the 0.3- 
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Fig. 3. Coadded EPIC spectra of the nucleus on 2002 Jun 7 
(blue) and Dec 8 (red), fitted with a constant two-temperature 
vmekal model plus a varying power-law component. 




channel energy (keV) 

Fig. 4. Coadded EPIC spectra of the ULX on 2002 Jun 7 and 
Dec 8, fitted with a disk-blackbody and a power-law, respec- 
tively. 



Firstly, we analysed the individual pn and MOS spectra from 
each observation. After ascertaining that they were consistent 
with each other, we coadded them with the method described 
in Page et al. (2003), in order to increase the signal-to-noise 
ratio. Thus, we obtained one combined EPIC spectrum of the 
nucleus and one of the off-nuclear source from each epoch. 
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12keVbandis * SxlO'^^ergs ' in June and 1.1 xlO'*^ ergs ' 
in December. The fitted temperatures of the optically-thin ther- 
mal components are kT^ = 0.43+° j° keV and kT2 = 0.98+^|° 
keV. Thermal plasma emission contributes for ^ 30-40% of 
the emitted luminosity in the 0.3-12 keV band and « 75- 
80% in 0.3-1 keV band. The power-law component has a slope 
T X 1.9, consistent both with an AGN spectrum and with typ- 
ical X-ray binary spectra. The emitted luminosity in the 2-10 
keV band is « 2.2 x lO"*" erg s^^ and =^ 3.8 x 10"*° erg s^^ in the 
two epochs: this is consistent with a current star-formation rate 
of a few Mo yr"! (Gilfanov et al. 2004). 

As an aside, we also searched for X-ray emission from 
NGC 7715, the tidally-interacting companion to NGC 77 14. 
Current star-formation from this galaxy is known to be neg- 
ligible, probably because of a relative lack of gas (Struck & 
Smith 2003). We did not significantly detect any source from 
each observation; however, when we combine both the June 
and December datasets, we detect a faint, possibly extended 
source at ^ 3-cr level in the EPIC image (Fig. 1). Its loca- 
tion (R.A. (2000) = 23'' 36"* 23*.6, Dec. (2000) = -1-02° 09' 32") 
puts it in the eastern gas/stellar countertail of NGC 7715 rather 
than in its nucleus. From the detected count rate, we estimate 
an average luminosity « 10^^ erg s'^ 

3.2. Off-nuclear ULX 

We then examined the coadded EPIC spectra of the off-nuclear 
source from the two observations. It is immediately clear that 
this source is fundamentally different from the starburst nu- 
cleus: its X-ray spectrum is featureless, with no significant evi- 
dence of optically-thin plasma emission. It is also clear that the 
source varied significantly between the two epochs. 

The source cannot be well fitted with a single spectral 
model in both observations. Limiting our choice to simple 
models, we find that a multicolor disk-blackbody with kTi^ = 
1.0 + 0.1 keV provides the best fit to the 2002 June spectrum 
(xl - 1.11; see Table 3 and Fig. 4). A simple power-law fit 
gives r = 2.5 + 0.2 (xl - 1-20), but is slightly improved by 
the addition of a blackbody component at kT = 0.66^^11 keV 
(xl - 1.15). On the contrary, a simple power-law model with 
F = 2.1;|;q j is an excellent fit to the 2002 December spectrum 
(xl = 0.86); adding a thermal component does not improve 
the fit. A multicolor disk-blackbody model is clearly rejected 
(xl = 1.29) for the second observation. We also tried fitting 
the two spectra with the Comptonisation model bmc (Shrader 
& Titarchuk 1999) (Table 3), confirming that the source was 
softer in June. 

Although any estimates of the emitted luminosity are very 
model-dependent, it is clear that the source has become brighter 
by a factor of ~ 2 between the two epochs, reaching a luminos- 
ity « 6 x lO"'" erg s"^ in 2002 December.^ Such variability on a 



^ The fact that both the ULX and the nuclear source showed an 
increase in luminosity from June to December prompted us to double- 
check that this effect is not due to systematic errors. The change in 
flux is at least an order of magnitude greater than any possible recip- 
rocal contamination between the two sources. Moreover, we verified 
that the luminosity and spectral variations are consistently seen in the 



Table 2. Best-fit parameters for the combined EPIC spectra of 
the nuclear source on 2002 Jun 7 and Dec 8. We assumed a 
Galactic column density ne.Gai - 4.9 x 10^*^ cm"^. We fitted the 
two spectra simultaneously, assuming that the thermal plasma 
component did not vary between the two epochs, but leaving 
the power-law component free to vary. The quoted errors are 
the 90% confidence limit. 



Parameter 


Value in Jun 02 


Value in Dec 02 


model: wabsoai x (wabsvm x (vmekali -i- vmekal2) + wabspo x po) 


nH,vm (10^' cm-2) 








«H.po (10^' cm-2) 








A:r™i (keV) 




43+" '" 

"■^•'-0.16 




fcr,„2 (keV) 




o.98:^t^ 




K.mi (10-5) 




3.6!?J 




Kyna (10-5) 




2.0!;:? 




r 


2.00!°f 




1 Ojr+0.19 
i.O->_0.19 


(10-5) 






6.0!1:^ 


C 




(0.8!°:^) ^ 




N 




(o.8!»:^)" 













Ne 




1 2*^-^ 

^•^-1.2 




Na 




(1.0) 




Mg 








Al 




(1.0) 




Si 




3 6*5.* 
-'•"-1.7 




S 




(1.0) 




Ar 




(1.0) 




Ca 




(1.0) 




Fe 




(1.0) 




Ni 




d.O) 




xl 




0.77 (122.4/160) 


fo.i-n (10-" ergcm-2 s^ 


-') ^Kl 




4 o+"-^ 


/0.3-1 (10-" ergcm-^ s" 


') o.8!°:i 




"•'-0.4 


L03-12 (10* erg s-') 


7.8 




10.5 


io.3-i2,vm(10*«ergs-i) 




2.9 




io.3-i2,po (10*« erg s-i) 


4.9 




7.6 


Lo.3-1 (lO-*" erg s-') 


3.5 




3.9 


L2-10 (10""* erg s-') 


2.2 




3.8 



We imposed the same abundance for these three elements. 



individual pn and MOS spectra; this rules out a systematic error in our 
coadding of the datasets. In the nucleus, we can clearly attribute the 
change to the power-law component only, while the thermal plasma 
emission is unchanged: this also confirms that it is not an instrumental 
artifact. As a further check, we extracted and analysed the spectra of 
the nearby QSO [HB89] 2333 + 019 (Fig. 1) for the two epochs. We 
find that the best-fit parameters (power-law with F = 1.65*q[J, and 
F = 1.79*Q Jq respectively) are consistent with an unchanged spec- 
trum. In both observations, the detected QSO flux in the 0.3-12 keV 
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2002 Jun 07 



5000 10' 
Time from start of observation (s) 



2002 Dec 08 

^ t \ L 

1.5x10^ 



Fig. 5. EPIC lightcurve for the ULX on 2002 Jun 7 and Dec 8, 
binned to 500 s. The background- subtracted count rate is the 
combined MOS plus pn count rate within an 11" extraction 
region; the encircled energy within that radius is ~ 59% for 
MOSl, and ^ 54% for MOS2 and pn. 

6-month timescale suggests that the size of the emitting region 
is ^ 0.2 pc. 

In addition to the long-term variability of the ULX flux, we 
looked for short-term variations during the two observations. 
We find that the background-subtracted 0.2-12 keV count rate 
is consistent with a constant level in the low state. However, 
variability is detected in the high state (Fig. 5), significant at the 
90% level according to the Kolmogorov-Smirnov test. Fitting 
the data with a sin curve gives a period of 6920;^^^"" s, and a 
semi-amplitude of 18%. However, longer observations are 
needed to ascertain whether the lightcurve has a real period- 
icity. No significant variability is detected from the nucleus 
within each of the two observations. 



4. Discussion 

4.1. Intermediate-mass BH or beamed emission? 

It was suggested (Papaderos & Fricke 1998) that the off-nuclear 
X-ray source could be emission from thermal plasma in a hot 
spot, due either to the collision of a fast outflow with the denser 
gas in the galactic bridge; or, to the infall of high-velocity 
clouds along the bridge onto the outer HI disk of NGC 7714. 
However, both those hypotheses are now ruled out by our 
XMM-Newton study. From the X-ray spectral properties in the 
two observations, and from the variability between the two 
epochs, we conclude that the off'-nuclear source is instead an 
accreting compact object. As such, with an extrapolated bolo- 
metric luminosity 1.5 x 10^' erg s"' in the high state, it is one 
of the brightest ULXs ever found. If the emission is isotropic, 
satisfying the Eddington limit would require a mass of « a few 
xIO^Mq for the accreting BH. 

The observed transition between a low/soft and a high/hard 
spectral state is opposite to what is generally observed in 

band was = 2.7x 10"'^ erg cm"^ s"' . We conclude that the variations 
seen in the ULX and the nucleus between 2002 June and December 
are real. 



Table 3. Best-fit parameters for the combined EPIC spectra of 
the ULX on 2002 Jun 7 and Dec 8. The quoted errors are the 
90% confidence limit. We assumed a Galactic column density 

"H.Gai = 4.9 X lO^" cm-l 



Parameter 



Value in Jun 02 Value in Dec 02 



model: wabsoai x wabs x (bb + po) 



nn (10^' cm-2) 


2 2+2-' 


1 5+°" 

^•-'-0.4 


kT^b (keV) 


66+" 

"•""-0.24 




^^bb (10-') 


4 4+4.2 
^- -3.5 




r 


2.6!«;^ 


2 l+o-^ 






6 r" 

O-i_0.9 



Xy 

/0.3-12 (10-'-' . 

Lo.3-,2(10*°ergs-') 



1.15(43.7/38) 
1 1+"' 
4.3 



0.86 (64.8/75) 

9 ^+0.2 
^■■-'-0.2 

6.6 



model: wabsoai x wabs x diskbb 



Mh (10^' cm-2) 
kT,, (keV) 
^dbb (10-3) 



07+" '" 

"■"'-0.07 

97+" 12 

"•^'-0.13 

6 2+" ' 

"•^-2.1 



< 0.12 

1 07+°'' 

^•"'-0.10 
7 4+2.6 
'■^-2.5 



Xy 



1-13 , 
,40 , 



io.3-i2(10''"ergs-i) 



1.11(44.4/40) 
1 0+" ' 
2.3 



1.29 (96.7/75) 
1.8!: 
3.8 



model: wabsoai x wabs x bmc 



n„ (10^' cm-2) 


<0.37 


34+- '' 


fcTbh (keV) 


34+" "* 

"■-^ -0.04 


22+"" 

"• -0.22 


r 


9 Q+0.4 
^■^-0.3 


1 9+"' 

l-5'_04 


log(A) 


> 0.10 


> 0.21 


^bmc (10-') 


Q 4+0 V 


14 5+^ " 



" erg cm ^ ' '^ 

i40 



io.3-i2(10''"ergs-') 



1.16(44.1/38) 

i.ii?' 

2.4 



0.87 (63.6/73) 
9 4+0.1 

5.1 



Galactic stellar-mass BH candidates, but similar to what was 
found for flie ULX Holmberg II X-1 (Dewangan et al. 2004) 
and for some ULXs in the Antennae (Fabbiano et al. 2003b). 
In the high state, the X-ray spectrum can be modelled by a pure 
power law with photon index F ^ 2; the conventional explana- 
tion for this component is inverse-Compton scattering of soft 
photons in a hot corona (Sunyaev & Titarchuk 1980). 

A thermal component (blackbody or disk-blackbody) is in- 
stead significantly detected in the low/soft state: it contributes 
for « (0.5-2.3) X lO""' erg s-' depending on the choice of spec- 
tral model; the short duration of out XMM-Newton observations 
does not allow us to put stronger constraints on the thermal 
emission. When the low-state spectrum is fitted with a disk- 
blackbody model, the color-temperature parameter kTi„ « 1 
keV (Table 3). If this parameter is identified with the effec- 
tive temperature at the inner boundary of a Shakura-Sunyaev 
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disk (Shakura & Sunyaev 1973), its high value is inconsistent 
with the BH mass inferred from the X-ray luminosity. High 
disk temperatures inconsistent with the standard disk model 
have been found in many other ULXs (Makishima et al. 2000). 
They have been explained (Mineshige et al. 2000; Watarai et 
al. 2001) with a different structure of the accretion flow at high 
accretion rates ("shm-disk"). Using the slim-disk model, our 
fitted value of Ti^ corresponds to a BH mass » 200Mo (in the 
Schwarzschild case), with a mass accretion rate a: 3 x 10"^ 
yr Another suggestion put forward to explain the high color 
temperatures in ULXs is a correction to the disk-blackbody 
spectrum via a hardening factor T\alT^s ~ 2.6 (Shrader & 
Titarchuk 2003). Finally, a standard accretion disk around a 
fast-rotating Kerr BH would extend closer to the event horizon 
and therefore would also have a higher inner temperature than 
predicted in the non-rotating (Schwarzschild) case (Zhang et 
al. 1997). However, this efl'ect would be significant only at high 
inclination angles (edge-on systems) and is unlikely to be the 
correct explanation for the high-temperature ULX disk spectra 
(Ebisawa et al. 2003). 

An alternative to the intermediate-mass BH scenario is the 
possibility that we are looking into a very narrowly collimated 
beam of emission (beaming factor ^ 50) from a stellar-mass 
X-ray binary (BH mass ^ 2OM0). Although the latter scenario 
cannot be ruled out a priori, it would be a very peculiar coinci- 
dence to find that the only strongly collimated source aligned — 
by chance — with our line of sight is located precisely at the 
junction between the coUisional stellar ring and the tidal bridge. 
A more general, quantitative argument against simple geomet- 
rical beaming is discussed in Davis & Mushotzky (2004), based 
on the consideration that the scattering region (thick disk or 
torus) responsible for the beaming can never be a perfect mir- 
ror; a large fraction of the radiation would be absorbed and 
re-radiated isotropically in the optical/lR. Hence, one would 
expect to find many bright (absolute magnitudes up to » -11) 
optical/IR point sources for every ULX. Such bright sources 
have never been seen in nearby galaxies. 

Alternatively, the observed ultraluminous flux might be due 
to relativistic beaming (Urry & Shafer 1984; Georganopoulos 
et al. 2002; Kording et al. 2002). In this scenario, the power- 
law component in the X-ray spectrum is the result of inverse- 
Compton scattering of soft photons by the electrons in a 
jet, moving with velocity f5 and a bulk motion factor F = 
(1 - ffy^^^. Relativistic Doppler boosting implies that the 
observed luminosity in a given band is enhanced by a fac- 
tor [F(l -/3cosff)Y'^. Assuming that the accreting source in 
NGC 7714 is a stellar-mass BH, and that the mechanical lumi- 
nosity in the jet is ^ O.lLEdd (eg. Fender 2001), Lorentz factors 
F ^ 3 and beaming angles ^ 5° are required to explain the 
observed luminosities in the high/hard state. However, this sce- 
nario would not explain the soft thermal component seen in the 
low state, which is unlikely to be beamed. 

In conclusion, the location, luminosity and spectral prop- 
erties of the ULX favour a massive (M > IOOMq) accreting 
object. Further constraints can come from time-variabiUty stud- 
ies. For example, a break in the power-density spectrum (PDS) 
at a frequency » (M/M©) Hz has been detected in some AGN 
and Galactic BH binaries (for recent examples see Markowitz 



et al. 2003; Uttley et al. 2002; Czerny et al. 2001), and from a 
ULX in NGC 4559 (Cropper et al. 2004). A preliminary PDS 
analysis of the ULX in NGC 7714 in its high state suggests fea- 
tures at ~ 1.5 X 10 Hz (also apparent from the lightcurve in 
Fig. 5) and ~ 2 x 10"'' Hz; hence, longer observations could be 
interesting. 

4.2. ULX formation in colliding galaxies 

One of the few statistically significant properties of the ULX 
population in nearby galaxies is that these sources are more fre- 
quently found in merging or tidally interacting galaxies (Swartz 
et al. 2003). It is not yet clear why. A simple suggestion is that 
the ULX formation rate is proportional to the star-formation 
rate (eg, Gilfanov et al. 2003), which is known to be higher in 
interacting galaxies. More specifically, ULX formation could 
be related to the clustered star-formation rate. 

The ULX in NGC 7714 is located on the outer stellar ring, 
at the junction with the tidal bridge. ULXs are often found in 
coUisional rings and tails of tidally disrupted systems, for ex- 
ample in the outer stellar ring of the Cartwheel galaxy (Gao et 
al. 2003). However, the two situations are fundamentally dif- 
ferent. In the Cartwheel, the outer ring is defined by an ex- 
panding wave of star formation (Higdon 1995) associated to 
a strong density wave in the gas, triggered by the coUisional 
perturbation; hence, ULXs in the Cartwheel ring are clearly as- 
sociated with active star formation and young stellar clusters. 
Here, instead, the off-centre collision between the two galax- 
ies produced a weaker dynamical perturbation (Struck & Smith 
2003), and the coUisional outer ring is only an expanding den- 
sity wave of old stars (Bushouse & Werner 1990), without on- 
going star formation: no Ha emission is observed along the 
ring (Gonzalez-Delgado et al. 1995; Smith et al. 1997). Radio 
21 -cm VLA observations also do not show any HI counterpart 
to the ring (Smith et al. 1997). 

Moderately active star formation and HI gas (total HI mass 
= 1 .6 X IO^Mq) are instead found along the tidal bridge (Smith 
et al. 1997). A string of young star clusters along the bridge is 
also visible in HST/WFPC2 images (Fig. 2; see also Struck & 
Smith 2003), although no obvious counterpart is seen in the 
ULX error circle (this is partly because a WFPC2 chip gap 
runs exactly across the circle). From the level of Ha emission 
in the maps of Smith et al. (1997) and Gonzalez-Delgado et 
al. (1995), the star-formation rate per unit area in the western 
part of the bridge (where it joins the stellar ring, near the ULX 
location) seems to be ~ 3-3.5 orders of magnitude lower than 
in the nuclear starburst region (total Ha luminosity from the 
bridge ^2x10^^ erg cm"^: Smith & Struck 2001). Assuming 
that the X-ray luminosity scales with the Ha emission and the 
star-formation rate (Ranalli et al. 2003; Buat et al. 2002) the 
extended X-ray emission component from the bridge would 
be ^ 10-'** erg cm"^, at least a factor of 10 fainter than our 
XMM-NewtonfEPlC detection limit. Despite the comparatively 
weaker star-forming activity, the bright ULX is located there 
rather than in the nuclear region. So, the effect of tidal interac- 
tions on ULX formation must be more complex than a simple 
increase of the star-formation rate. 



R. Soria and C. Motch: A variable ultra-luminous X-ray source in the colliding galaxy NGC 7714 



7 



We can clearly rule out the possibility that the ULX be a 
massive runaway binary formed in the nucleus and then mi- 
grated to the outer ring. If, as suggested by its luminosity, the 
total mass of the binary system is > 25Mq, we expect kick ve- 
locities < 10 km s~^ (Zezas & Fabbiano 2002). At a distance 
of « 4 kpc from the nucleus, it would have taken the ULX 
> 4 X 10** yr to get there, more than the lifetime of any mas- 
sive donor stars and possibly more than the age of the nuclear 
starburst itself (Lanfon et al. 2001). 

Despite the high column density of atomic hydrogen 
10^1 cm-^: Smith et al. 1997), CO (1-0) observations failed 
to detect the tidal bridge (Smith & Struck 2001). This can be 
interpreted as evidence for the low metal abundance of the gas. 
The galactic nucleus is also known to be metal-poor (~ 0.2- 
O.4Z0: Garcia- Vargas et al. 1997). Other bright ULXs have 
been found in interacting galaxy systems, located at or near 
large, high-density HI or H2 complexes, but with only a weak 
level of local star formation. For example, the ULX M8 1 X-9 
is associated (Wang 2002) with a molecular cloud or "proto- 
galaxy" near Hohnberg IX, with a gas mass ~ IO^Mq and star- 
formation rate ^ 10 "^M® yr ' (Brouillet et al. 1992; Henkel 
et al. 1993). In fact, the last close encounter between M81 
and M 82 seems to have left a number of tidal tails and de- 
bris (Yun et al. 1994), stripped from the metal-poor outer part 
of M 81, which are now evolving into tidal dwarf galaxies, with 
metal abundances ~ 0.1-0.4Zo (Boyce et al. 2001; Makarova 
et al. 2002). 

As pointed out by Pakull & Mirioni (2002), stellar evo- 
lution in metal-poor environments may favor the formation 
of ULXs, because the mass-loss rate in the radiatively-driven 
wind from the massive O-star progenitor is much reduced 
(M„ ~ Z^-^^: Vink et al. 2001; see also Bouret et al. 2003). 
This leads to a more massive stellar core, which may then col- 
lapse into a more massive BH. Hence, it may explain the for- 
mation of BHs with masses up to « SOM© and isotropic lumi- 
nosities X 10"*" erg s"'. Tidal dwarfs, bridges and tails stripped 
from the outer disks of colliding galaxies may provide an en- 
vironment for the formation of molecular clouds and stars out 
of metal-poor gas. However, normal stellar evolution processes 
alone would be unable to explain accreting BHs with isotropic 
luminosities « lO"*' erg s"^ such as the one seen in NGC 7714. 

BHs with masses of ~ lO^-lO^M© may be formed in the 
core of young (age ^3x10* yr), massive star clusters, due to 
the Spitzer instability, runaway core collapse and merger of the 
O stars (Portegies Zwart & McMillan 2002; Rasio et al. 2003). 
Deeper optical observations of the ULX field in this galaxy, and 
a more precise determination of its position^, are needed to in- 
vestigate the possible association with massive young clusters. 

Alternatively, BH remnants in this mass range are thought 
to be left over by primordial-metallicity (Population III) stars 
(Madau & Rees 2001). Most of these BHs could reside in 
galactic halos: assuming halo masses of 6 x lO'^Mg and 2 x 
lO^^'Mo for NGC 7714/15 respectively (Struck & Smith 2003), 
and using the results of Islam et al. (2003), one can estimate the 
presence of ~ 200 such intermediate-mass BHs in this system. 
In an undisturbed galactic halo, accretion from the low-density 

^ The system has never been observed by Chandra. 



interstellar medium would make them too faint to be detected. 
Could galactic collisions and mergers create the conditions for 
some of these BHs to become bright (up to their Eddington 
luminosity ~ 10"*' erg s '), by providing a fuel supply? This 
could in principle occur either via tidal capture of a donor star, 
or via accretion from a dense molecular cloud. The tidal cap- 
ture timescale can be comparable with the dynamical timescale 
of the NGC 7714/15 interaction only for stellar densities ^ 10^ 
stars pc"^ (Fabian et al. 1975; Zezas & Fabbiano 2002), much 
larger than in the collisional stellar ring. Hence, this process 
should be ruled out. 

On the other hand, let us consider Bondi-Hoyle accretion 
from an intergalactic molecular cloud. Assuming an efficiency 
77 = 0.1, an isotropic luminosity « lO"*' erg s"^ requires an 
accretion rate M ^ 10^' g s"' « 2 x IO^^Mq yr"'. The Bondi- 
Hoyle accretion rate (Bondi & Hoyle 1944; Mirabel etal. 1991) 
is 

M« 2.4x 10"(M/Mo)^(n/cm^)(v/10kms-^)-^ gs'^ 

where M is the mass of the accreting BH, n is the number den- 
sity of the molecular gas at large distances from the BH, and 
V is the velocity of the compact object relative to the cloud. 
Choosing M - 10^ Mq and v = 10 km s \ H2 number densi- 
ties « 2 x 10"* cm"^ are required to explain the observed lumi- 
nosity. Such densities are typical of molecular cloud cores (eg. 
Plume et al. 1997), with sizes ~ 0.1 pc. The total molecular 
gas mass in the NGC 7714 bridge is likely to be ~ IO^^-IO'^Mq, 
depending on the metal abundance (Smith et al. 1987; Smith & 
Struck 2001). Hence, the filling factor of dense molecular cloud 
cores near the western end of the tidal bridge is only ~ 10 
The probability that a primordial BH could encounter one such 
cloud appears to be low. Addressing these issues in more de- 
tails is beyond the scope of this work. 

4.3. Nuclear activity 

Our XMM-Newton study shows that the X-ray spectrum of 
the starburst nucleus contains a variable power-law compo- 
nent in addition to non-variable thermal plasma emission 
(Table 2). This suggests that there is probably a single ac- 
creting point-source in the nuclear region contributing a sig- 
nificant fraction of the power-law emission, at least 3 x 10'*" 
erg s '. It had been speculated in the past that NGC 7714 
might contain a low-luminosity AGN, but it was then shown 
(O'Halloran et al. 2000) that a normal starburst can account 
for the IR/optical/UV properties. However, these new obser- 
vations suggest the presence either of an AGN, or of another 
ULX. Other cases of galaxies with variable AGN activity un- 
detected from their optical colors or spectra are discussed in 
Davis & Mushotzky (2004) and references therein. A more de- 
tailed study of the nuclear starburst, combining the new X-ray 
data with the multiwavelength results of Langon et al. (2001), 
and Gonzalez-Delgado et al. (1995, 1999), is left to further 
work. 

Finally, we note that the NGC 7714/15 system is one of 
the examples of alignments between a galaxy and two quasars 
discussed by H. Arp and collaborators (Stocke & Arp 1978). 
Within that scenario, it has been suggested (Arp et al. 2004) 
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that ULXs and quasars might be related to the same physical 
process of ejection from galactic nuclei. We point out for com- 
pleteness that the ULX in NGC 77 14 is also aligned with the 
nucleus and the two quasars. However, we have no elements to 
speculate that this represents more than a coincidence. 

5. Conclusions 

We have used XMM-Newton to study the interacting galaxy 
system NGC 77 14/15. We have reported here the main prop- 
erties of the two brightest sources: the starburst nucleus and an 
off-nuclear ULX. The X-ray spectrum of the off-nuclear source 
suggests that it is an accreting BH, and rules out the possibility 
that it is due to thermal-plasma emission from a hot spot, as pre- 
viously speculated. Its X-ray flux varies by a factor of 2 over a 
six months' interval; the source appears softer in the low state, 
unlike the typical behaviour of Galactic BH candidates but in 
agreement with the behaviour of many ULXs. Its spectrum in 
the low/soft state can be fitted by a disk-blackbody model with 
kTin ~ 1 keV: this is inconsistent with a Shakura-Sunyaev disk, 
but can be explained with a slim-disk model. 

In the high state, its emitted isotropic luminosity is » 6 x 
lO'**' erg s"' in the 0.3-12 keV band, implying a bolometric lu- 
minosity a! 1.5 X 10^' erg s"' from a reasonable extrapolation 
of the power-law spectrum (photon index F x 2). BH masses of 
~ a few 1Q^-1Q^Mq would be required to satisfy the Eddington 
limit. Furthermore, variability on timescales of a; 2 hr is de- 
tected in the high state. 

The ULX is located at the junction of the tidal bridge (con- 
sisting of gas and young stars) with the collisional outer ring 
(consisting of an old stellar population, with no gas). We have 
pointed out that ULXs are often found in tidally interacting 
systems, associated with metal-poor molecular clouds, tidal 
dwarfs, or HI structures formed in the galactic collision. 

The nucleus has an X-ray luminosity ^ lO"*' erg s"' in 
the 0.3-12 keV band. Thermal plasma emission contributes for 
» 3 X 10"*° erg s"', constant over the two observations, and 
is probably extended (marginally resolved in the EPIC/MOS 
images). A point-like power-law component contributes for 
~ 5x lO'*" erg s"' and a; SxlO"**^ erg s"' in the two observations. 
The power-law component in the X-ray spectra of starburst nu- 
clei is generally due to unresolved high-mass X-ray binaries. 
The amount of variability in our case implies that one single 
source contributes for at least 3 x 10'*° erg s"'. This suggests 
that there is either a hidden AGN or another ULX in the nu- 
clear region. 

Acknowledgements. We thank Manfred Pakull for the preparation of 
the observations and for discussions. 
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